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ABSTRACT 

Pulsar efficiency, defined as the ratio of the pulsed bolometric luminosity to the spin-down 
power, is calculated to be ~ 15% (averaged over the spin-dipole inclination angle, ranging between 
« 50% for the aligned and w 10% for the orthogonal). We also estimate the characteristic photon 
energy and argue that our results agree with the Fermi pulsar catalog - in a sense. 

1. Introduction 

Pulsars create magnetospheres around them- 
selves (Goldreich and Julian 1969, Ruderman and 
Sutherland 1975). Full quantitative description of 
a magnetosphere requires a kinetics input - we 
need to know how the charges and photons are 
produced. 

It turns out that the main overall characteris- 

o 

tics of the pulsar are calculable without any de- 
tailed knowledge of the kinetics. Here we calcu- 
late the pulsar spin-down power L s d and the pul- 
sar efficiency e, defined as the ratio of the pulsed 
bolometric luminosity to L sc [- 

Our results, valid to some 10% accuracy, are as 
follows. The efficiency is 
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where 9 is the spin-dipole angle. The spin-down 
luminosity is 
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where /i is the magnetic dipole moment and f2 is 
the angular velocity. 

2. Observations 

(i) The Fermi pulsar catalog (Abdo et al 2010) 
lists 46 pulsars. Only 10 are known to have e < 
0.1. These might be out of the beam. 

(ii) Our theory (§3) in the present form should 
not be used for light curves and spectra - kinet- 
ics seem to matter. But a crude estimate of the 



Fig. 1. — The damping region of a pulsar with 
the spin-dipole angle 8 = 45°. Ohm's law P = 
\/ p 2 + 10(E 2 + B 2 )/r 2 , r is the spherical radius. 
The star rotates around the z axis at half the speed 
of light. Shown are the isolines of E 2 — B 2 (where 
positive) in the plane y = 0. The snapshot is 
taken near the time when the maximum of E 2 — B 2 
passes through the plane y = 0. The resolution is 
100 3 . 
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characteristic photon energy is possible (curvature 
emission in an electric field ~ magnetic field at the 
light cylinder, radius of curvature ~ light cylin- 
der radius): E ~ dhe-^L^Qi ~ SL^PmsGeV, 
where L34 is the spin-down power in units of 
10 34 erg/s and P ms is the pulsar period in ms. This 
formula is OK (to factor 3) for all the catalog pul- 
sars but 14. For these 14 pulsars, the formula gives 
a frequency some 3-10 times higher than observed. 
All these 14 pulsars either have or may have small 
efficiency, and may be out of the beam. 

3. Electrodynamics of Massless Charges 

We calculate the pulsar using the so-called Elec- 
trodynamics of Massless Charges (EMC, Gruzinov 
2012). EMC describes the motion of light charges 
in strong electromagnetic fields. Due to strong 
radiation damping, the dynamics becomes "Aris- 
totelian" , with the charge velocity rather than ac- 
celeration given by the field values. Positive and 
negative charges move at the speed light in the 
directions 
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Here the scalar E and the pseudoscalar B are 
the proper electric and magnetic fields defined by 

Bl - El =B 2 - E 2 , B E = B • E, E > 0. (4) 

To calculate the electromagnetic field, we nu- 
merically solve Maxwell equations. To solve 
Maxwell equations we need to know the current. 
Let p± be the absolute values of the charge den- 
sities of positive and negative charges. Then the 
EMC equation (J3J gives 
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where p = p+ — p- and P = p + + p_ . Since p = 
V-E, we almost have an Ohm's law (an expression 
for j in terms of E and B): 
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Equation ([6]) is not quite an Ohm's law, because 
P is unknown (we only know that P > \p\). 

We have checked that when the neutral plasma 
density P — \p\ becomes large enough, in practice 



for P — \p\ > \dB\, the simulation results become 
insensitive to the actual value of P. This happens 
in the following way. In most of the volume, the 
E and B fields are force-free, pE + j x B = 0. 
There is also a small volume where the fields are 
not force- free. In this volume some of the Poynting 
flux emitted by the star gets dissipated. For larger 
P, the non-force-free volume shrinks, but the to- 
tal dissipated power remains fixed. The damping 
region is shown in Fig. 1. 

4. Appendix: Numerics 

(i) We use uniform cubic grid in a cube. E com- 
ponents are at the corresponding (parallel) edges 
of the grid, B components are at the correspond- 
ing (orthogonal) faces of the grid. We use outgoing 
boundary conditions at the faces of the cube. We 
solve the evolutionary problem 



-V x E, <9 t E = V x B - j. 
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(ii) We have a rotating star at the center of the 
cube. This means that inside some sphere the cur- 
rent j is given by the Ohm's law in a rotating frame 
plus an external rotating current density. The ex- 
ternal current density is axisymmetric, with the 
symmetry axis of the current inclined with respect 
to the rotation axis. We get the dipole radiation. 

(iii) Then we add the Ohm's law dU) outside the 
star. At this stage we have to regularize equations 
([7|) by adding small diffusivities (77AB and rjAE 
with small 77) . The Poynting flux emanating from 
the star is the spin-down power L sc i- The Poynting 
flux leaving the simulation cube is (1 — e)L s d- 
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